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A new simplified approach based on equilibrium modeling is proposed in this work to describe the cor¬ 
relations among syngas species experimentally observed in a pilot scale downdraft biomass gasifier oper¬ 
ated with different feedstocks (biomass pellets and vine prunings). The modeling approach is based on 
experimental evidence on the presence of devolatilization products in the syngas and fluctuations of syn¬ 
gas composition during stationary operation, accounted for by introducing two empirical parameters, a 
by-pass index and a permeability index. The simplified model correctly reproduces the correlations 
among the main syngas species (including methane and ethylene) resulting from experimental data of 
pilot tests with different feedstocks and under a wide range of operating conditions. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The prediction of the composition of the syngas produced by 
biomass gasifiers is an important part in the development of mod¬ 
eling tools, which can support the understanding of the gasifier 
phenomenology as well as the assessment of the gasification 
performance. 

This task requires an accurate representation of several phe¬ 
nomena involved in a gasification system, such as devolatilization, 
heterogeneous and homogeneous reactions, heat and mass trans¬ 
fer. In most cases, the development of complex models may result 
inapplicable, due to the lack of information and data needed to run 
the simulations. 

As far as fixed bed gasifiers are concerned, a simple approach to 
compute the syngas composition relies on equilibrium models (Zai- 
nal et al., 2002; Altafini et al„ 2003; Jarungthammachote and Dutta, 
2007; Baluand Chung, 2012; Barman et ah, 2012). These models of¬ 
ten consider a stoichiometric combustion of the biomass in the oxi¬ 
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dation zone followed by an equilibrium stage. As observed by Balu 
and Chung (2012) these models often underestimate the methane 
content in the syngas and to this purpose Barman et ah (2012) 
introduced a deviation from the equilibrium in the methane forma¬ 
tion reaction to fit its prediction with the experimental values. 

Equilibrium models cannot account for the presence of ethylene 
in the syngas. Although its concentration is usually low compared 
to other major species, ethylene is usually more abundant than 
other hydrocarbons (Wander et ah, 2004; Simone et ah, 2012a) 
and it is considered an indicator of biomass devolatilization (Boro- 
son et ah, 1989; Fagbemi et ah, 2001). 

Model predictions are usually compared with average syngas 
composition, which is clearly satisfactory when the model is used 
to evaluate gasification performance. However, experimental tests 
at pilot scale (Simone et ah, 2012a,b) show that the syngas compo¬ 
sition is characterized by a dynamic evolution, affected by operat¬ 
ing parameters and bed properties. 

To be reliable, a simplified gasification model should take into 
account at least the variation of the syngas composition and be 
able to describe the presence of the most important hydrocarbons 
in the syngas. 
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Nomenclature 



OCi 

devolatilization yield of the ith component (-) 

Pn 

pressure at the nozzle outlet (mm H 2 0) 

b 

by-pass index (-) 

A Pj 

pressure difference between the atmosphere and the 

Cij 

stoichiometric combustion coefficient of the ith compo- 


annular jacket of the gasifier (mm H 2 0) 


nent in reaction; (-) 

A P n 

pressure difference between the atmosphere and the 

Cpi 

specific heat of ith component (J moG 1 K 1 ) 


nozzle outlet of the gasifier (mm H 2 0) 

s 

heat loss coefficient (-) 

PM, 

molecular weight of ith component (g moG 1 ) 

Sij 

stoichiometric gasification coefficient of the ith compo- 

Pcj 

combustion reaction rate (kmol hr 1 ) 


nent in reaction; (-) 

Pgi 

gasification reaction rate (kmol h ') 

Ah cj 

enthalpy of reaction c; 

rij 

tar cracking and hydrocarbon reforming reaction rate 

m bio 

biomass loading rate (kg h ') 


(kmol nr 3 h ') 

nf em 

ith component molar flux after drying and devolatiliza- 

tii 

stoichiometric tar cracking and hydrocarbon reforming 


tion (kmol tr 1 ) 


coefficient of the ith component in reaction tj (-) 

n oxi 

ith component molar flux from the oxidation zone 

■j’devo 

devolatilization zone temperature (K) 


(kmol tr 1 ) 

■j'Oxi 

oxidation zone temperature (K) 

nj ed 

ith component molar flux from the reduction zone 

yby-pass 

by-pass volume in the oxidation zone (m 3 ) 


(kmol tr 1 ) 

^moi 

biomass moisture content (-) 

n by-pass 

ith component molar flux from the by-pass (kmol h 1 ) 

X 

permeability index (-) 

Po 

atmospheric pressure (mm H 2 0) 



Pj 

pressure in the annular jacket (mm H 2 0) 




In this work a new simplified approach based on equilibrium 
modeling is proposed which reproduces the variation of the syngas 
composition observed in a pilot scale downdraft gasifier operated 
with different feedstock. 


2. Methods 

2.1. Experimental setup 

The experimental tests are carried out with a pilot scale down- 
draft biomass gasifier (described in detail in Simone et al. (2012a)). 
The plant is operated slightly below atmospheric conditions due to 
a fan-blower positioned at the end of the gas clean-up line, which 
drives air to enter the gasifier through four nozzles positioned in 
the throated section of the gasifier. This oxygen rich section is 
called oxidation zone, since combustion reactions occur there. 
The section above the oxidation zone is oxygen-free but in the 
proximity of the oxidation zone the temperature is sufficiently 
high for biomass drying and devolatilization to take place. Hetero¬ 
geneous gasification reactions occur in the section below the 
throat, called reduction zone. 

The produced syngas moves upward from the bottom of the 
gasifier in an external annular jacket and enters the clean-up sys¬ 
tem where the syngas is scrubbed with water and filtered before 
being analyzed. 

The gasification plant is equipped with two on-line instruments 
for gas sampling and analysis: a micro-gas chromatograph (micro- 
GC) and a Fourier Transform Infra-Red Spectrometer (FTIR). 

The micro-GC (Agilent 3000) allows analyzing hydrogen, light 
hydrocarbons and permanent gases in around 4 min. The micro-GC 
is equipped with two independent channels based on an injector, a 
column and a thermal-conductivity detector (TCD). The first channel 
is based on a Molsieve 5A column, using argon as mobile phase and is 
suitable for the separation of hydrogen, oxygen, nitrogen, methane 
and carbon monoxide. The second channel is based on a PLOT U col¬ 
umn using helium as mobile phase and suitable for the separation 
of carbon dioxide, ethane, ethylene and acetylene. The micro-GC is 
equipped with a pump for gas sampling and a membrane filter. 

FTIR measurements are performed using a Bruker Tensor 37 
FTIR spectrometer, installed in a specifically developed pressurized 


cabinet. The apparatus is equipped with a pump for gas sampling 
including valves and manometers for pressure control and a heated 
transfer line. A heated gas cell (stainless steel body, 10 cm optical 
pathlength, 25 ml volume, 200 °C maximum temperature) is used 
for gas analysis. 

Both the analyzers are connected to the plant after the syngas 
clean-up, consequently the gas at the sampling point is free of 
water vapor and the compositions are referred to the diy syngas. 

The gasifier behavior during operation is monitored using a per¬ 
meability index, defined as: 

„ AP, 

X "A T, 

where: 

• APj is the pressure difference between the atmospheric pres¬ 
sure (P 0 ) and the pressure in the annular jacket of the gasifier 

(PjY- 

A Pj = P 0 - Pj 

• AP n is the pressure difference between the atmospheric pres¬ 
sure (P 0 ) and the pressure at the nozzle outlet inside the gasifier 

(P n ): 

A P„ = Po ~ Pn 

The permeability index displays a remarkable influence on the 
gasifier operation, affecting the syngas productivity and operative 
stability (Simone et al., 2012a). In particular, air adduction and syn¬ 
gas flow-rate increases with the bed permeability. 

2.2. Descriptive approach 

Several experimental tests with the pilot scale downdraft bio¬ 
mass gasifier were performed using different types of biomass as 
feedstock (Simone et al., 2012a,b). In particular, seven tests are 
considered for the analysis performed in this paper. Wood pellets 
were used in Test 1, Test 2 and Test 3; a mix of wood pellets and 
sunflower cake pellets was used in Test 4 and Test 5; vine prunings 
were used in Test 6 and Test 7. The syngas was assayed by micro- 
GC in Tests 1-7, while the FTIR was used in Test 7 only. 

The analysis of the experimental results obtained in different 
tests highlighted the fluctuation of the main chemical species in 
the syngas (CO, C0 2 , H 2 , CH 4 , C 2 H 4 ). For instance, Fig. la shows 
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Fig. 1. Example of microGC (a) and FTIR (b) signals obtained during Test 7. 


an example of the microGC data obtained during Test 7 in the gas¬ 
ification of vine prunings. In the same test, the fluctuations are 
even more evident from the FTIR signals reported in Fig. lb, due 
to the higher sampling and detection frequency. As can be ob¬ 
served by comparing data in Fig. la and b, the assays carried out 
with both methods (microGC and FTIR) provide similar measure¬ 
ments for all relevant gas species (with the exception of hydrogen, 
not detected by FTIR), with a relative difference of around 6% (aver¬ 
age of all data reported in Fig. 1 ). The composition variability, typ¬ 
ical of any test carried out in the pilot scale gasifier, is a 
characteristic feature of the gasification process, independent of 
the feedstock, syngas sampling frequency, and the assay method 
(microGC or FTIR). 

In spite of the general stability in terms of operating conditions, 
these results show a continuous dynamics of the gasifier. As dis¬ 
cussed below, the analysis of experimental data shows that, de¬ 
spite their apparent randomness, these fluctuations can be 
associated to rather clear correlations between syngas species. This 
dynamics may be partly ascribed to the variation of the permeabil¬ 


ity index X, which was found to vary throughout the experimental 
tests (Simone et al., 2012a). Changes in the bed permeability affect 
the air adduction to the biomass bed, which in turn affects the 
instantaneous equivalence ratio and the gasifier temperature pro¬ 
file, and therefore the syngas composition. 

The experimental data confirm also that the ethylene concen¬ 
tration in the syngas is not negligible (Simone et al., 2012a,b), thus 
suggesting a direct influence of devolatilization products on the fi¬ 
nal syngas composition. Ethylene and methane are products of pri¬ 
mary and secondary pyrolysis reactions (Boroson et al., 1989; 
Fagbemi et al., 2001). This step influences the final syngas compo¬ 
sition when the volatiles produced from pyrolysis are not com¬ 
pletely oxidized in the oxidation zone of the gasifier. This may be 
due the occurrence of preferential channels within the bed which 
allow pyrolysis gases to by-pass the oxidation and reduction zone 
(Tarhan, 2003; Yamazaki et al., 2005) and may result in inefficien¬ 
cies of the oxidation zone in abating tars. 

These experimental evidences (i.e. composition correlations and 
devolatilization products in the syngas) are not accounted for by 
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Devolatilized biomass: Char, Tar, H 2 , CO, C0 2 , CH 4 , C 2 H 4 , H 2 0 
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Fig. 2. Pictorial representation of the modeling approach. 


available equilibrium models (Zainal et al., 2002; Altafini et al., 
2003; Jarungthammachote and Dutta, 2007; Balu and Chung, 
2012; Barman et al., 2012), and a new approach is proposed in this 
work to overcome their limitations. This approach relies on two 
parameters: 

- The aforementioned permeability index (X), governing the air 
adduction to the gasifier. 

- The by-pass index (b), accounting for the contribution of devol¬ 
atilization to the final syngas composition. 


The approach is pictorially represented in Fig. 2. The gasifier is 
divided into different connected sections. 

In the upper section of the gasifier the biomass is completely 
dried and devolatilized into char, tar, water and permanent gases. 
Devolatilization is schematized as: 


BIOMASS —> cqFb + a 2 CO + a 3 CH 4 + a 4 C 2 H 4 -f a 5 C0 2 + a 6 H 2 0 
-|- a 7 Tar+ otgChar 


The components in the devolatilization reaction (gas species: 
hydrogen, carbon monoxide, methane, ethylene and carbon diox¬ 
ide; water; tar; char) are those characterizing any process in the 
gasifier. They are referred to in the model with a numerical index 
(i) ranging from 1 to 8, corresponding to their stoichiometric coef¬ 
ficient in the devolatilization reaction. 

Char is considered as pure carbon. The biomass ultimate analy¬ 
sis used in the calculations is that of wood pellets reported in Si¬ 
mone et al. (2012a). The coefficients a, are assigned on the basis 
of the work of Fagbemi et al. (2001), the tar raw formula is then 
calculated on the basis of the elemental balances (C, H, O). 

The mass balances for the gas species (i = 1-5), tar (i = 7) and 
char (i = 8) in the upper section read as follows: 


a,- 

PMi 


ttlbioO ^moi) 


0) 


The mass balance for water (i = 6) must include the feedstock 
moisture content: 


,devo 

'h 2 


^— m bi0 [a 6 (l Xrnoi) “F ^rnoij 

™ H 2 0 


( 2 ) 


The lower section of the downdraft gasifier is schematized with 
two parallel reactors: 


- A mixed reactor representing the ideal behavior of the gasifier. 


- A by pass section, accounting for the fraction of devolatilization 
products which by-pass the oxidation and reduction zone. 


The by pass index is defined as the fraction of devolatilization 
products entering the by-pass section. 

The devolatilization products in the mixed reactor are burned 
with air in the oxidation zone. This process is represented as a se¬ 
quence of combustion reactions which occur according to the oxy¬ 
gen availability. Oxidation of the pyrolysis products is allowed to 
consume the available oxygen in the following sequence of 
descending order of reactions rates: hydrogen (reaction cl), ethyl¬ 
ene (reaction c2), methane (reaction c3), carbon monoxide (reac¬ 
tion c4). If more oxygen remains after oxidation of the pyrolysis 
gas products, char (reaction c5) starts to be oxidized. The products 
of hydrocarbon combustion are water and carbon monoxide which 
is further oxidized to carbon dioxide in a following step. A similar 
approach was used by Sharma (2011 ). The char conversion through 
heterogeneous oxidation decreases with increasing flux of volatiles 
in the oxidation zone. 

In the oxidation zone, tar is assumed to be completely con¬ 
verted through cracking reactions (reaction c6) to give H 2 , CO, 
CH 4 , C 2 H 4 and char. 

The mass balances for the gas compounds (i = 1-5), water (i = 6) 
and char (i = 8) in the oxidation zone read as follows: 

n° x ‘ = (1 - b) ^ n ? evo + J^CjjRcj'j ( 3 ) 

where (j = 1-5) are the rates (kmol h ’) at which the gas species 
and char are oxidized depending on the availability of oxygen, k c6 is 
the rate at which tar is cracked, and c,, are the stoichiometric coef¬ 
ficients for the reaction of the ith component in the jth reaction, 
negative for reactants and positive for products. 

The oxidation products (water, carbon dioxide, carbon monox¬ 
ide, nitrogen and residual char and volatiles in excess) enter the 
reduction zone, represented as a chemical equilibrium reactor 
where the following reactions occur: 

C + H 2 O^H 2 + CO (gl) 

C + C0 2 ^2C0 (g2) 

C0 + H 2 0^H 2 + C0 2 (g3) 

The mass balances for the gas compounds (i = 1-5), water (i = 6) 
and char (i = 8) in the reduction zone read as follows: 

n[ ed = n° xi + ( 4 ) 

j=i 


where the (j = 1-3) are the rates at which the gas species, water 
and char are consumed and/or produced when the reactions gl -g3 
are at equilibrium, and gy are the stoichiometric coefficients for the 
reaction of the ith component in the jth reaction, negative for reac¬ 
tants and positive for products. 

The air adduction to the reactor is a function of the permeability 
index (X) which is varied according to the typical ranges observed 
in the experimental tests. The air flow-rate to the gasifier, and thus 
the equivalence ratio, increases with increasing permeability 
index. 

The temperature of the oxidation zone is computed on the basis 
of the heat balance over the throated section of the gasifier: 

5>° xi c pi (r xi - r d ™) = (1 - 5)J2 R cAh cj (5) 

1=1 j=l 


where S is a lumped parameter, estimated using the experimental 
temperature profile in the gasifier, accounting for heat dispersion 
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from the oxidation zone to the other sections of the reactors (dev¬ 
olatilization and reduction zones) and heat losses through the reac¬ 
tor walls. The inlet temperature to the oxidation zone (T devo ) is 
estimated on the basis of the experimental observations on the 
temperature profile in the reactor. 

The temperature of the reduction zone (and thus of the equilib¬ 
rium reactor) is related to the oxidation zone temperature on the 
basis of the experimental observations of Simone et al. (2013). 

The other parallel reactor takes into account the portion of dev¬ 
olatilization products which by-pass the oxidation and reduction 
zone. This is a completely stirred tank reactor where tar cracking 
and reforming reactions can occur. 

Tar cracking (reaction tl) is represented according to the work 
of Boroson et al. (1989), whose yields are used to calculate the 
amount of carbon monoxide, carbon dioxide, methane, ethylene 
and hydrogen generated by the cracking reaction. The kinetics used 
in the model is that of Liden et al. (1988). Hydrocarbon reforming 
(reaction t2) was treated as in Tinaut et al. (2008). 

The mass balances for the gas compounds (i = 1-5), water (i = 6) 
and tar (i = 7) in the by-pass zone read as follows: 

^by-pass = j n devo + (joxi) yby-pass 

V J=1 

where r tJ are the reaction rates (krnol m~ 3 hr 1 ) of tar cracking 0 = 1) 
and hydrocarbon reforming 0 = 2), and ty are the stoichiometric 
coefficients of species i in reaction j. 

The temperature of the reactor is that computed for the oxida¬ 
tion zone and the residence time is calculated according to the con¬ 
tinuity equation of the gas phase assuming incompressible flow. 

Finally the streams from the two parallel reactors are mixed to 
represent the syngas produced by the gasifier. 

The by-pass index (ranging from 0 to 1) is defined as the frac¬ 
tion of devolatilization products entering the by-pass section. High 
values of by-pass index correspond to high levels of the contribu¬ 
tion of devolatilization to the final syngas composition, represent¬ 
ing extensive channeling. The ideal equilibrium condition in the 
gasifier is approached when the by-pass becomes negligible. 

Calculations were carried out with a fixed biomass loading rate 
of 50 kg/h, equal to the nominal throughput of the gasifier and 
varying independently the permeability index (X) from 2 to 12 
(as a typical range of variation in pilot tests, Simone et al., 
2012a) and the by-pass index ( b ) from 0 to 0.9. 

Model equations are a set of non linear algebraic equations 
which was solved using the family-of-lines method implemented 
by the gPROMS platform (PSE, London). 

3. Results and discussion 

3.3. Experimental results 

Although the syngas composition appears to randomly fluctuate 
during stationary operation of the gasifier (see Fig. 1), clear trends 
and correlations among chemical species can be identified. For in¬ 
stance, Fig. 3 shows some composition plots for the gasification of 
pelletized biomass obtained by microGC. 

An inverse proportionality between the carbon dioxide and car¬ 
bon monoxide volume fractions can be observed in Fig. 3, which 
can be related to both variation of the oxygen availability and tem¬ 
perature of the reduction zone. The oxidation yield of char and gas¬ 
eous species to carbon dioxide and the temperature in the gasifier 
increase with increasing oxygen content. Higher temperatures re¬ 
sult in higher carbon monoxide production yield according to the 
Boudouard and water gas shift equilibria. 

The hydrogen volume fraction is somewhat independent of the 
carbon monoxide or carbon dioxide concentrations. The bulk of 


experimental volume fractions data ranges from 15% to 20% but 
this variation cannot be clearly related to variation of the other 
gaseous species. 

Clear trends can be observed by plotting the methane and eth¬ 
ylene concentrations against the carbon dioxide volume fractions. 
The two graphs are similar, showing exponential relationships be¬ 
tween the two hydrocarbons and carbon dioxide. Opposite trends 
are observed by plotting the methane and ethylene volume frac¬ 
tions against the carbon monoxide volume fractions (not reported). 
This highlights also an important relationship between the meth¬ 
ane and ethylene content. Notably there is a direct proportionality 
between the two species though the data in the plot appear quite 
scattered. 

Fig. 4 shows the composition plots obtained in two gasification 
tests of vine prunings (Tests 6 and 7). In Test 7 data were obtained 
with the FTIR analyzer which allowed a higher sampling frequency 
(7 s), with respect to the microGC (4 min). 

Data from Tests 6 and 7 slightly differ for the carbon monoxide- 
carbon dioxide ranges. Test 7 is characterized by a lower content of 
carbon dioxide, whereas the two sets of data show very similar 
relations between methane and ethylene. 

The plots reported in Fig. 4 are similar to those observed in the 
gasification of biomass pellets, but in this case the trends are 
clearer and less scattered. In addition, the ranges of variation are 
more limited. For instance, in the gasification of vine prunings 
the methane and ethylene contents range from 1.6% to 3.3% and 
from 0.2% to 0.6%, respectively. In the gasification of biomass pel¬ 
lets the two species range from 1.5% to 4% and from 0.2% to 1%. 

With the exception of different CO-C0 2 range (which can be as¬ 
cribed to different test history), the comparison of the FTIR and mi- 
cro-GC data does not show remarkable differences in the CO-C0 2 
relationship and the C 2 H 4 -CH 4 relationship, confirming that the 
trends are independent of the sampling frequency. On the contrary, 
the C 2 H 4 -C0 2 and the CH 4 -C0 2 plots show more scattered plot 
when the FTIR is used. This may suggest that the CO-C0 2 and 
C 2 H 4 -CH 4 relationships are stronger than those between oxidized 
species and hydrocarbon species. 

Given the limited scattering of the data obtained with vine pru¬ 
nings, it is possible to estimate average values for the ratios CO¬ 
CO;. (1.722 ±0.410 in Test 6 and 1.776 ±0.420 in Test 7) and 
C 2 H 4 -CH 4 (0.166 ±0.012 in Test 6 and 0.168 ±0.020 in Test 7). 
Noteworthy the C 2 H 4 -CH 4 ratios observed in two different tests 
are almost identical and characterized by a very low standard devi¬ 
ation, thus indicating that the formation of the two species is clo¬ 
sely linked. 

The variation in the syngas composition observed in the two 
sets of tests may be partly ascribed to the variation of bed perme¬ 
ability, which affects the instantaneous equivalence ratio and the 
gasifier temperature profile. From this point of view, vine prunings 
allow for gasification conditions which are more stable than bio¬ 
mass pellets due to the different mechanical properties of the feed¬ 
stock. The higher stability of the gasification conditions of vine 
prunings may explain the less scattered plot and the clearer trend 
obtained for this feedstock. 

However, the sole variation of bed permeability and thus of air 
flow to the gasifier cannot explain the relationships between the 
gaseous species in the syngas. 

In particular, the close relationship between methane and eth¬ 
ylene in the syngas cannot be explained according to the kinetic 
and equilibrium models present in the literature. 

In gasifier models (Di Blasi, 2000; Tinaut et al., 2008) meth¬ 
ane is usually produced from biomass pyrolysis, tar cracking 
and carbon gasification with hydrogen. Ethylene is usually not 
explicitly considered but included in lumped species as light 
hydrocarbons (Thunman et al., 2001) and generated from bio¬ 
mass pyrolysis and tar cracking. The traditional representation 
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Fig. 3. Comparison of model results and experimental data for the gasification of pelletized biomass (Tests 1-5). 


of the downdraft gasifier phenomenology considers that all the 
volatiles produced in the pyrolysis zone burn in the oxidation 
zone. The presence of methane is then explained according to 
the gasification of char with hydrogen, but ethylene is usually 
not taken into account. 

However, given the straight relationship between ethylene and 
methane observed in the experimental tests, the two species are 
likely to be produced from the same process. As a consequence, 
the importance of the carbon gasification with hydrogen to pro¬ 
duce methane seems limited (as pointed out for atmospheric con¬ 
ditions by Di Blasi, 2000). 

3.2. Calculation results 

3.2.1. Parametric study 

The descriptive approach presented before was used to inter¬ 
pret the experimental relationships among the syngas species. 

The calculated trends are remarkably influenced by the calcu¬ 
lated average temperature of the reduction zone. Fig. 5a shows 
the calculated reduction zone temperature plotted against the 
by-pass index for different values of the permeability index. 

It can be seen that the average temperature of the reduction 
zone increases with increasing permeability index, due to the high¬ 
er equivalence ratio attained. 


The effect of the by-pass index is not as linear. At X = 2 an in¬ 
crease in the by-pass index leads initially to higher reduction tem¬ 
perature, which reaches a maximum at b = 0.4 and then starts to 
decrease. At higher X the maximum shifts towards lower values 
of b. At X = 4 the maximum is at b = 0.2. This is due to representa¬ 
tion adopted for the oxidation zone (see Fig. 2). 

At X = 2 and b = 0, almost all the carbon monoxide produced by 
devolatilization and hydrocarbon oxidation is unconverted due to 
lack of oxygen. Therefore the temperature is low since much of 
the energy retained in the biomass is released in this step. The in¬ 
crease of the by-pass index reduces the amount of hydrocarbons to 
the combustion section, therefore there is more oxygen available 
for carbon monoxide combustion which leads to a temperature in¬ 
crease in the oxidation and, consequently, reduction zone. When 
carbon monoxide is totally converted (b = 0.4 when X = 2) the max¬ 
imum temperature is reached. A further increase of the by-pass in¬ 
dex leads to reduction of the carbon monoxide availability, with 
oxygen starting to consume char. As a result the temperatures de¬ 
creases, since less heat is released per mole of oxygen consumed. 
The value of b at which carbon monoxide is depleted and the max¬ 
imum temperature is attained decreases with increasing X. 

In Fig. 5a the average temperatures of the reduction zone range 
between 743 and 988 K. Not all these temperatures are likely to oc¬ 
cur in a biomass gasifier. Temperatures below 873 K seem too low 
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0 0.2 0.4 0.6 0.8 1 


By-pass index 

Fig. 5. Calculated (a) reduction zone temperatures and (b) tar content plotted 
against the by-pass index for different permeability index. 


if compared with temperature profile of the gasifier (reported in Si¬ 
mone et al. (2013)), though this does not cover the whole reduction 
zone. Therefore very high values of the by-pass index seem unli¬ 
kely to occur. 

The same considerations can be drawn for the oxidation zone 
temperature. The calculated temperatures for the oxidation zone 
range from 991 to 1317 K. Again the values obtained with high 
by-pass index are outside the experimental measures reported in 
Simone et al. (2013), which were between 1173 and 1473 K. 


3.2.2. Comparison of experimental and model data 

Figs. 3 and 6 compare the results of the descriptive approach 
with the composition plots obtained from the experimental tests 
on biomass pellets and vine prunings, respectively. The calculated 
results are reported as different series for different values of by¬ 
pass index. Each series is composed by a set of data which repre¬ 
sent the variation of the permeability index. Series with by-pass in¬ 
dex 0 and 0.2 show only three data corresponding to permeability 
index of 2, 3 and 4. Higher values of X cannot be computed for 
these low by-pass index configurations since, at higher equivalence 
ratio, the equilibrium conditions require a char consumption high¬ 
er than that supplied to the gasifier. 

When the by-pass index is set to zero the gasifier is represented 
as an ideal equilibrium reactor. This representation is far from the 
experimental evidence, since no methane or ethylene is present in 
the syngas, irrespectively of the permeability index and thus of the 
equivalence ratio, and the calculated hydrogen values (25%) are 
largely overestimated both for biomass pellets and vine prunings. 

When the by-pass index increases to 0.2, methane and ethylene 
are produced in the syngas. The hydrogen and carbon dioxide con¬ 
tents decrease, while the carbon monoxide content increases. 

This is a consequence of the higher temperature of the reduc¬ 
tion zone attained at b = 0.2, which favors the carbon monoxide 
production at the expenses of carbon dioxide and hydrogen, 
according to the water gas shift reaction. 
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Fig. 6. Comparison of model results and experimental data for the gasification of vine prunings (Test 6). 


Increasing b to 0.5 leads to higher methane and ethylene con¬ 
tents. At the same time there is a further reduction in the hydrogen 
content, which at this point is mostly due to the dilution of the 
hydrogen rich gas stream produced in the ideal section with the 
hydrogen lean produced in the devolatilization. It can also be seen 
that the effect of the permeability index on the composition be¬ 
comes less important, as at this point the final syngas composition 
is mostly independent of this parameter. 

Starting from this point, further increases of the b index show a 
clear trend. The final syngas composition becomes more similar to 
the devolatilization gas containing more ethylene and methane. 
The reduction of the temperature of the reduction zone, reduces 
the amount of hydrogen and carbon monoxide in favor of carbon 
dioxide produced from char oxidation. In fact, at very high values 
of b, the ideal section acts as a char combustor producing almost 
only C0 2 , which cannot be reconverted to CO due to the low tem¬ 
peratures of the reduction zone. 

The comparison of the calculated compositions with the exper¬ 
imental data is good for vine prunings (Fig. 6) starting from b = 0.5 
(X = 2-4) to b = 0.6 (whole X range), except for a slight underesti¬ 
mation of the methane content. Notably compositions with b rang¬ 
ing between 0.5 and 0.6 are those which fit all the composition 
plots for all the values of permeability index. Outside this range 


not all the volume fractions of the different gaseous compounds 
are properly computed. 

For vine pruning tests, the independent variation of the two 
parameters seems to provide a good representation of the compo¬ 
sition variation, owing to the rather close gathering of experimen¬ 
tal data. 

The same range of b index seems to provide a satisfactorily rep¬ 
resentation of the composition plots obtained for biomass pellets 
gasification. In this case the data are more scattered than with vine 
prunings, so the use of the two parameters cannot fully represent 
the gasifier behavior, especially for the variation of the ethylene 
and methane contents. However, the model seems to provide a 
good description of the variation of the syngas composition even 
with a larger set of data. In particular it provides a reasonable 
explanation for the exponential trend exhibited by the C 2 H 4 -C0 2 
and CH 4 -C0 2 plots and the quasi-independency of hydrogen on 
the permeability index. 

The amount of devolatilization gases by-passing the bed in¬ 
creases, while the temperature of the reduction zone decreases, 
with increasing by-pass index. In addition, a lower temperature 
of the reduction zone implies a lower conversion of char into gas. 
Both these effects results in higher fractions of methane and ethyl¬ 
ene in the produced syngas. 
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Table 1 

Proximate and ultimate analysis of different solid residues produced by the gasifier. 
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Test 

VM (% dry) 

FC(% dry) 

Ash (% dry) 

C (% daf) 

H (% daf) 

N (% daf) 

Test 3 

7.65 

83.70 

8.65 

84.79 

0.44 

0.64 

Test 5 

5.96 

86.14 

7.90 

84.60 

0.20 

0.65 

Test 6 

13.1 

78.5 

8.4 

83.80 

0.90 

1.00 

Test 7 

9.3 

81.0 

9.7 

82.10 

0.90 

1.00 


The poor dependency of hydrogen on the permeability index is 
related to the opposite variation of the heterogeneous gasification 
with water and the water gas shift reaction, which are improved at 
higher and lower temperatures, respectively. The two reactions 
vary with the temperature in such a way that the increase of the 
temperature of the reduction zone slightly affects the net hydrogen 
production. 

Since most of the experimental data are fitted with values of b 
in the range 0.5-0.6, it can be concluded that the gasifier never 
operates according to the ideal representation and that there is al¬ 
ways a certain extent of devolatilization gas contributing to the fi¬ 
nal composition. 

This has an important implication on the tar content in the syn¬ 
gas. Since tars are meant to be destroyed in the throated section, 
characterized by high temperatures, the occurrence of channeling 
allows for the production of tar which cannot be abated in the gas¬ 
ifier. Therefore, tar content in the syngas increases with increasing 
by-pass extent. This is clearly highlighted by the model. Fig. 5b 
shows the tar content in the syngas plotted against the by-pass in¬ 
dex for different permeability index. 

It can be seen that increasing the by-pass index the tar content 
increases dramatically. In the range b = 0.5-0.6 the tar content is 
still below 0.01%, but a further increase of the by-pass index 
greatly enhances it due to the combined effect of the higher by¬ 
passed fraction and the lower temperature which reduce the effect 
of the cracking reactions. 

Notably the tar content calculated by the model is consistent 
with the typical values reported by Stevens (2001). Taking the 
values of Stevens (2001) as reference, the gasifier seems to operate 
already in the lower range when the gasifier is operated at b = 0.6 
and to produce a rather low tar content when b is 0.5. Therefore if 
compared to the literature values the tar production of the down- 
draft gasifier seems acceptable. Nevertheless the gasification unit 
has to be equipped with a suitable clean-up system in order to 
comply with the end-user requirements reported in Hasler and 
Nussbaumer (1999). In order to optimize the performance of the 
downdraft gasifier it is necessary to relate the by-pass index to 
the operating parameters of the gasifier. Yamazaki et al. (2005) 
have correlated the tar yield from a downdraft gasifier with 
superficial gas velocity, showing that at high gas superficial 
velocities the tar yield increases due to channeling and short resi¬ 
dence time. 

Downdraft gasifier are generally operated at the highest syngas 
flow-rate to maximize the power generation. However, this implies 
high gas superficial velocities in the downdraft gasifier and thus a 
certain extent of by-pass. Instead, good operation practices should 
be developed to attain the highest syngas flow-rate provided that 
by-pass phenomena are still below a certain level. These proce¬ 
dures should be defined on the basis of the end-user requirements, 
the clean-up capabilities and the gasifier configuration. 

In this work, the presence of devolatilization products in the 
syngas is related to channeling. Experimental works indicated that 
the solid residue produced by fixed bed gasifiers usually exhibits 
very high carbon contents (Klinghoffer et al., 2012; Ma et al., 
2012) and in some cases volatile matter (see Table 1). These results 
confirm the non ideal behavior of downdraft gasifiers. The volatile 
matter may be tarry compounds condensed on the cooler char bed 
of the reduction zone, thus confirming the occurrence of channel¬ 


ing. On the other hand, biomass devolatilization may be still taking 
place in the reduction zone. The discharge frequency of the bottom 
grid is often increased to cope with high pressure drops through 
the bed. This practice reduces the residence time of the biomass 
in the gasifier, and may lead to non-completely devolatilized bio¬ 
mass to reach the reduction zone. According to this hypothesis 
the gasifier performance should not be evaluated only in terms 
of optimal pressure drop, as this practice may lead to high tar con¬ 
tents as well as lower biomass conversions. 

The degree of by-pass and the bed permeability can be con¬ 
trolled to some extent by satisfying the gasifier specifications for 
feedstock size distribution, density and morphology. Fine particles 
reduce the void fraction and the pore size in the gasifier bed, thus 
reducing the bed permeability and leading to lower syngas produc¬ 
tivity. On the other hand, oversized particles (and elongated parti¬ 
cles in particular) may lead to bridging and uneven distribution of 
the feedstock in the gasifier, causing by-pass of the oxidation and 
reduction zones for a portion of the gas flow in the reactor. The 
results presented in Figs. 3 and 6 do not show a clear difference 
in the range of by-pass index fitting the experimental data for 
biomass pellets (Fig. 3) and vine prunings (Fig. 6). More research 
work is required to relate the feedstock properties (morphology, 
size distribution and density) to the permeability index and by¬ 
pass index, introduced in this work as parameters characterizing 
the performance of downdraft biomass gasifiers. 

4. Conclusions 

The extended equilibrium model presented in this paper allows 
the prediction of syngas composition, including devolatilization 
products (methane and ethylene), in downdraft biomass gasifica¬ 
tion. The simulations, validated using pilot-scale experimental 
data, suggest that degree of by-pass and bed permeability are 
key macroscopic factors determining the syngas composition and 
therefore the gasifier performance. More research work is required 
to assess the effects of operating conditions and feedstock proper¬ 
ties (morphology, size distribution, density) on these factors. The 
by-pass index may depend on gasifier configuration and size, and 
validation should then be extended to other case studies to im¬ 
prove the model prediction capability. 
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